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Abstract. The computations of numerical relativity make use of (3 + 
1)— decompositions of Einstein field equations. We examine the concep- 
tual characteristics of this method; instances of compact-star binaries 
are considered. The preeminent role of the geodesic motions is empha- 
sized. 



PACS 04.20 - General relativity. 
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1. Introduction. Insect. 2 we recall properties and geodesies of the Gauss- 
ian-normal coordinate system. In sect. 3 we emphasize that the physical 
results derived from any solution of Einstein field equations must be inde- 
pendent of the adopted reference frame. In sects. 4 and 5 we give a resume 
of the properties of the (3 + l)-decompositions of Einstein equations, and we 
point out that: i) the results of the approximate computations obtained with 
the employment of the gravitational energy-pseudotensor have an illusive 
value; ii) coordinate-system independence of the results is not proved; Hi) 
no mathematical theorem of existence supports the approximate computa- 
tions of numerical relativity. Sect. 5bis: the particles of a discrete "cloud of 
dust" describe geodesic lines - and therefore no GW is emitted by them: an 
emblematic instance. Sects. 6 and 7 contain some comments on the numeri- 
cal computations concerning the binaries composed of two "Schwarschildian 
mass-points" and of two neutron stars; an analysis is given of the notions 
of ADM-mass. Appendix A: The Einstein equations in a Gaussian-normal 
frame. Appendix B: The equations of the standard (3 + l)-decomposition. 
Appendix C: The four-dimensional world as "a mass of plasticine" (Weyl), 
and some interesting consequences. 



2. - The concept of Gaussian-normal coordinate system is due to Hilbert 
PQ. Landau and Lifshitz called it synchronous [SJ. Let us consider in the 
four-dimensional world a three-dimensional space S3 such that every line in 
it is space-like; (x 1 , x 2 , x 3 ) be the point-coordinates in 53. Starting from 
each point (x 1 , x 2 , x 3 ), we trace the geodesic lines which are orthogonal to 
53. They will become particular time- like lines if we report on them the 
time x° = ct as the proper time r. It follows easily that 
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(1) ds 2 = -(dx ) 2 +g a pdx a dxP , (q,/3 = 1,2,3) . 

By virtue of the very construction of the three-dimensional space x° = 0, 
the quadratic form g a p dx a dx 13 is necessarily a definite positive one. 

It is interesting to write the differential equations of all the geodesic lines 
of the metric ([1]). The Lagrangian £: 

(2) £ := -(x ) 2 + g a p(x°,x 1 ,x 2 ,x 3 )x a xP = Ac 2 , 

where the overdot denotes a derivative with respect to an affine parameter 
a, and A is a constant, gives the Lagrangian equations 

/o\ -o , 1 9g a p . a .R 

(3) x ^ = ; 



(3') 

g ai x« + ^x°x« + ^x«x»--^x"x>> = , (7,5 = 1,2,3). 



1 dg a/3 p 



£ is a first integral of eqs.(pH3j; A is negative,zero, positive for time- 
like, null, space-like geodesies, respectively. If a = r, we have A = — 1. 
The geodesies orthogonal to x° = are characterized by = x a = x a , and 
a = t. 

The solutions of eqs.fpHS" 7 ]) are commonly interpreted, if a = t, as the 
trajectories of test-particles. However, they represent also the geodesic tra- 
jectories of the material elements of a "cloud" , which generates the Gaussian 
field. 

There are infinite expressions for the Gaussian- normal ds 2 , which are ob- 
tained by one of them with: i) any Lorentz transformation of the spacetime 
coordinates, ii) any transformation of the space coordinates. 

Landau and Lifshitz [2j write the Einstein field equations in a Gaussian- 
normal frame; then, they separate space - and time - derivatives, thus 
obtaining a particular (3 + l)-formalism. All the operations of raising and 
lowering of the indices and the covariant derivatives are performed with 
respect to the space metric g a (3(x°, x 1 , x 2 , x 3 ). They put K a p := dg a p/dt, 
and in their equations (99,10), (99,11), (99,12) (see App. A) the symbol 
Pa denotes the Ricci tensor with respect to the i/iree-dimensional space. It 
is immediate to divide eqs. (99,10), (99,11), (99,12) in constraint equa- 
tions, containing g aj s and dg a /s/dt, and in evolution equations, containing 
d 2 g a p/dt 2 . It can be proved [3] that 



(4) 



d 2 g a p _ F al 3 



where g is the determinant of the Gaussian metric matrix, and F a s is an 
entire rational function of the g a ps, of the dg a p/dfs, of the second spatial 
and spatiotemporal derivatives of g a p. (Of course, F a p depends also on the 
matter tensor). 
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For the solution of the Cauchy problem it is necessary to give the values 
of g a p and dg a p/dt at all the space points (x 1 ,^ 2 ,^ 3 ) at an initial instant 
t = t - 



3. - In principle, the Cauchy problem could be solved for any physical sys- 
tem in any system of coordinates, thus giving origin to an implicit, particular 
(3 + reformulation of the problem. 

As it was emphasized by Hilbert pQ , as a rule any solution of Einstein field 
equations must satisfy the following criterion: all the physical results which 
we deduce from it must be independent of the coordinate system in which 
the solution is expressed. Of course, this does not exclude the existence of 
interesting physical phenomena - as, e.g., the gravitational redshift of the 
spectral lines -, whose entity depends on the chosen reference system. An 
emblematic instance of independence of the coordinate system is represented 
by the geodesic motions of the elements of a continuous "cloud of dust" . 



4. - The numerical relativity makes use of two (3 + 1 ^decompositions of 
the Einstein field equations [4a), b)}. The standard decomposition is (see, 
e.g., p.41 of [4a)]): 



(5) ds 2 = -a 2 dt 2 + 7ife (dx j + (3 j dt) (dx k + (3 k dt) , 

where (j, k) = (1, 2, 3), and c = G = 1. The matrix of the components of 
metric ([5]) is: 



(6) 



In Appendix B we report the field equations of the standard decomposi- 
tion. 

If a = 1 and the three-vector (3 3 is zero, eq. ([5]) gives a Gaussian-normal 
ds 2 . At p. 37 of [4a)] we read: "The freedom to choose these four gauge 
functions a, /3 , f3 2 , /3 3 completely arbitrary embodies the fourfold coordi- 
nate degrees of freedom inherent in general relativity ..." The lapse function 
a reflects the freedom to choose the sequence of the space-like hypersurfaces, 
and the shift vector f3 3 reflects the freedom to relabel the spatial coordinates 
on each of the above hypersurfaces in an arbitrary way. This means that 
the covariance of a (3 + l)-decomposition is restricted to the continuous 
transformations of of the space coordinates X j X j X , and to the freedom to 
choose the initiali-hypersurfaces. Consequently, the (3 + 1)- decompositions 
are not equivalent to the Einstein field equations, which possess the general 
covariance with respect to all the continuous transformations of the space- 
time coordinates 3 . The (3 + l)-decompositions are contrary to 
the spirit of GR. And in them a vast class of coordinate systems is not 
employable. 
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5. - In the applications to various problems of the formalisms of the (3 + 1)- 
decompositions of Einstein equations the authors choose suitable particular- 
izations of the metric (J5]H([6|). In these numerical computations an important 
role is commonly played by the gravitational energy-pseudotensor, which is 
a tensor only with respect to linear coordinate transformations, and there- 
fore its intervention renders illusive the physical value of any result. In- 
deed, it is not sufficient to remark that the total pseudoenergy and the total 
pseudomomentum of the system form a four-dimensional vector in the very 
distant Minkowskian spacetime. Further, the authors omit to prove that 
the obtained results are independent of the adopted reference system. No 
mathematical theorem of existence is given: the numerical solutions presup- 
pose that "suitable boundary conditions and initial data are chosen so that 
these solutions do indeed exist" ([4a)], p. 21). 

5bis. - Articles and treatises of numerical relativity and (3 + ^-decomposi- 
tions of the Einstein equations contain many considerations and compu- 
tations concerning the properties and the generation of GWs. We have 
given several demonstrations that GR, if properly understood, excludes the 
physical existence of GWs [3] - and the experience corroborates our proofs. 
Therefore, we do not discuss here the numerical-relativity methods regard- 
ing the GWs. We give instead the simple proof that also the motions of 
the particles of a discrete "cloud of dust" are geodesic. First of all, accord- 
ing to a method by Infeld [6] we write the Einstein equations, in a generic 
coordinate system x°, a; 1 , a; 2 , a; 3 , using tensor densities: 





^9tiiyV^gR = (ji, v = 0, 1, 2, 3) , (c = G = 1) . 

lf5(x—£) := ^(a; 1 — £ t 1 )5(x 2 — £, 2 )5(x 3 — £ 3 ) is the three-dimensional Dirac's 
distribution, and £ 1 ,£ 2 ,£ 3 are the coordinate of a particle, we put - if the 
"cloud" is composed of s particles: 



(8) 



p p 

(9) J=g T^=m(t) Hi_ 1L_ S ( z _ £) 



ds ds 

now, if the world-lines of the particles never intersect, it is not difficult to 
verify that the equations of motion of the particles: 

(10) 

P =i 

where the semicolon denotes a covariant derivative, are the differential 
equations of geodesic lines. Consequently, no GW is generated by our 
"dust"; a result very general, because there exists no limitation for the 
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values of the particle kinematical-elements (velocities, accelerations, time 
derivative of the accelerations, etc.). 

We see that the gravitational self-force theory, which is based on a wrong 
analogy with the electromagnetic self-force of a charge, does not make any 
physical sense. And the post- Newtonian approximations which concern a 
discrete "cloud of dust" cannot give GWs. 

6. - With regard to the binary stars composed of two Schwarzschildian 
mass-points (see, e.g., Chapt. 13 - "Binary Black Hole Evolution" - of 
[4a)]), we observe that nobody has succeeded in proving the theoretical ex- 
istence of two Schwarzschildian mass-points, and therefore the approximate 
solution of Einstein equations obtained by numerical computations, which 
would describe this system, has a very doubtful value. Further, we remark 
that the black-hole interpretation of the Schwarzschildian point-mass solu- 
tion of Einstein equations has a mathematically unfounded basis (see, e.g., 
sects. 3b, 3c of |5j). At any rate, the curvature ("hard") singularity at 
r = 0, which characterizes the gravitational field of a point-mass in a cur- 
rent consideration of the standard (Hilbert-Droste-Weyl) form of the metric, 
can be simply removed if one adopts the original Schwarzschild's coordinate 
system, or Brillouin's system, for which the metric has only a "soft" singu- 
larity at r = 0. Numerical stratagems, as the BH-excision, the moving 
puncture method, etc., are superfluous. 

7. - The notion of ADM-mass (ADM: Arnowitt-Descr-Misner) is cur- 
rently employed in the numerical computations [4a)]. Now, the definition 
of the ADM-mass is founded on the properties of the gravitational energy- 
pseudotensor (see sect. 5), and therefore the physical meaning of an ADM- 
mass is rather uncertain. 

Chapt. 16 of [4a)] regards the "Binary Neutron Star Evolution". The ap- 
proximate computations make use of the above mentioned notion of ADM- 
mass. We limit ourselves to emphasize that the time evaluations of the 
inspiral and merger phases of the considered binaries and of the evolutions 
of the amplitudes h+ and h x of hypothesized GWs (computed with the dis- 
putable quadrupole formula) do not possess a physical value. Indeed, these 
time intervals would make sense only if we could prove their independence 
of the adopted coordinate system. 



We conclude that the approximate computations of numerical relativity 
are essentially self-referential, because no experimental, or observational, 
proof and no theorem of mathematical existence supports them. Further, 
they do not even satisfy the criterion according to which the physical value 
of the results depends on their independence of the reference frame. 
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APPENDIX A 



The Einstein field equations written in a Gaussian-normal coordinate sys- 
tem (see sect. 2 and Landau and Lifshitz [2]) are (a,/3 = 1,2,3): 



(Al) 



_ J_ dK%_ 1 a 

K °~ 2c dt + 4 @ 



*E<L (T° - - T) 



(A2) 



(4 



„4 a 



(A3) 



rv l 



1 







2cy^g dt 



\T) ; 



here: K a /s := ® 9 £f ; P« is the three-dimensional Ricci tensor; the semicolon 
denotes a covariant derivative with respect to the three-dimensional metric 
g a/3: {a, (3 = 1,2,3). 



APPENDIX B 



We give the basic formulae of the standard (3 + l)-decomposition of the 
Einstein field equations [7J. 

Index notations: ...) = (1,2,3); (a, 6, ... = 0,1,2,3). The extrinsic 
curvature Kij = \/—g 00 Ti°j of a hypersurface measures how much normal 
vectors to the hypersurface differ at neighboring points; K = 7^ K %3 = 
Kj. The symbol Dj denotes the covariant derivative with respect to 7^. 
Normal vector: n a = (—a, 0,0,0); n a = (a -1 , — a~ l j3 l , — a~ l j3 2 , — a _1 /3 3 ). 
The matter tensor T ab is the fluidodynamical energy-tensor, with various 
polytropic equations of state (EOSs); the matter source terms are: 

g = n a n b T ab ; S i = - 7 « n a T aj ; Sij = lialjb T ab - S = 7^' S i3 . 

Metric (c = G = 1): 

(Bl) ds 2 = -a 2 dt 2 + 7 y (dx* + /3Mt) (dx j + /3 J dt) . 

Two constraint equations (R is the Ricci scalar): 

(B2) R + K 2 - Kij K ij = 16vrp ; 

(B3) Dj (K ij - f j K) = 8*5? . 

Evolution equation for the space metric 7^: 
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(B4) 



dt 



Evolution equation for the extrinsic curvature K, 



(B5) 
dK, 



it 



Of 



a {K^ - 2K ik Kj + KKij) ~ DiDja - 8vra 



Hj(S- q) 



dx k 



dx'i 



dx l 



APPENDIX C 



We read in Weyl [8] a penetrant analysis of a fundamental property of the 
general-relativistic coordinate systems. He wrote [46)]: "... the concept of 
relative motion of several bodies has, as the postulate of general relativity 
shows, no more foundation than the concept of absolute motion of a single 
body. Let us imagine the four-dimensional world as a mass of plasticine tra- 
versed by individual fibers, the world lines of the material particles. Except 
for the condition that no two world lines intersect, their pattern may be 
arbitrarily given. The plasticine can then be continuously deformed so that 
not only one but all fibers become vertical straight lines." It is clear that 
this consideration implies that there exists always a coordinate transforma- 
tion, which allows us to pass from any coordinate system for which some 
bodies are in motion to a co-moving coordinate system for which all these 
bodies are at rest. 

Now, no class of privileged coordinate systems exists in GR, and any 
physical effect must be frame independent [Ij. Consequently, the fact that 
bodies at rest cannot generate gravitational waves has a general significance: 
no coordinate system exists for which the motions of the bodies generate 
gravitational waves. 

Remark that these considerations hold for the general case in which both 
gravitational and non- gravitational forces are present. For a different proof, 
founded on the Einstein field equations, that all the general-relativistic mo- 
tions can be geodesically described, see our paper of ref. 

An immediate corollary: the gravitational field of a body, whose motion 
is geodesic, moves en bloc with the body, and is propagated instantaneously. 
In special relativity we have a partial analogue: the static-electromagnetic 
fields created by an electric charge in a rectilinear and uniform motion (a 
Minkowskian geodesic motion) move along with the charge, and are propa- 
gated instantaneously [10j . (This corresponds perfectly to the results found 
by an observer in a rectilinear and uniform motion, who travels with respect 
to the charge at rest). 
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This fact has been considered a paradox by some physicists, who have 
tried to get rid of it with the gratuitous surmise that it is hold only for 
infinite spatio-temporal motions of the electric charge. 

An instantaneous propagation of a field is not always in contradiction 
with the theory of relativity. - 
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